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Aims. The transiting extrasolar planet WASP-12 b was found to be one of the most intensely irradiated exoplanets. It is unexpectedly 
bloated and is losing mass that may accrete into the host star. Our aim was to refine the parameters of this intriguing system and search 
for signs of transit timing variations. 

Methods. We gathered high-precision light curves for two transits of WASP-12 b. Assuming various limb-darkening laws, we gen- 
erated best-fitting models and redetermined parameters of the system. Error estimates were derived by the prayer bead method and 
Monte Carlo simulations. 

Results. System parameters obtained by us are found to agree with previous studies within one sigma. Use of the non-linear 
limb-darkening laws results in the best-fitting models. With two new mid-transit times, the ephemeris was refined to BJDtdb = 
(2454508.97682 ± 0.00020) + (1.09142245 ± O.OOOOOOaS)^. Interestingly, indications of transit timing variation are detected at the 
level of 3.4 sigma. This signal can be induced by an additional planet in the system. Simplified numerical simulations shows that a 
perturber could be a terrestrial-type planet if both planets are in a low-order orbital resonance. However, we emphasise that further 
observations are needed to confirm variation and to constrain properties of the perturber. 
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1. Introduction 

WASP-12 b is one of about 100 extrasolar transiting planets 
known so far. The orbital plane of such a planet is aligned in 
such a way that the planet passes across the disk of its host star, 
temporarily blocking a fraction of stellar light. That effect is ob- 
served as a flux drop. Combining spectroscopic and photomet- 
ric data gives a unique opportunity to determine planetary mass 
and radius, hence density, a key parameter for studying internal 
structure of an exoplanet. 

Hebb et al. (2009 ) have found that WASP-12 has a transit- 
ing planetary-mass companion. The Jupiter-mass planet orbits 
the host star with a period of about 1.09 day causing transits 
with a depth of 14 milli-mag (mmag) and duration of 2.7 hours. 
The high-accuracy photometric follow-up observations of tran- 
sits were performed with the robotic 2-m Liverpool Telescope 
allowing the mid-transit time to be determined with a precision 
of ~20 s (Hebb et al. 120091 ). These observations allow the mass, 
the planet radius, and the orbital semi-major axis to be deter- 
mined as Mb = 1.41 ± 0.10 Ml, = 1.79 ± 0.09 R], and 
a - 0.0229 ± 0.0008 AU, respectively. The central star is found 
to be a late-F type dwarf whose age is estimated to be 2 + 1 Gyr 

The high effective temperature of the central star (Tetj - 
6300!^°° ^"'l ^hort orbital period make WASP-12 b one of 
the most intensely irradiated exoplanets. This results in a high 
equilibrium temperature of the planet equal to 25 16 + 36 K when 
assuming isotropic reradiation and a zero-Bond albedo. Hebb et 
al. (12009 ^ note that the internal source of energy may also be 
needed to explain the unexpectedly large radius of the planet. 



The orbit of such a close planet is expected to be circularised 
on short time scales, but surprisingly the orbital eccentricity of 
WASP-12 b was found to be nonzero (e = 0.049 + 0.015). Thus, 
dissipation of tidal energy has been pointed out as a mechanism 
heating up the planetary interior and bloating the planetary ra- 
dius (Miller et al. .2009 ). However, further observations showed 
that WASP-12 b's o rbit is highly circular with e = 0.017+° °[^ 
(Husnoo et al. '2010'). This result indicates that the eccentricity- 
driven tidal energy source is negligible. 



* Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA), operated jointly by the Max-Planck Institut 
fiir Astronomie and the Institute de Astrofisica de Andalucia (CSIC). 



Li et al. (120101) predict that the planet may be losing mass 
at a rate of ~10"^ Mj yr"' by exceeding its Roche lobe. The 
planetary gas is expected to fall towards the host star through 
Lagrangian point Li and form an optically thin accretion disk. 
Lai et al. (2010) find by modelling that an accretion stream can 
have a significant projected area, even comparable to the planet's 
projected area. The transfer of metals may enhance the apparent 
stellar metallicity. This effect should be easy to detect because 
WASP-12 is expected to have a very shallow convective zone. 
Using spectropolarimetric observations, Fossati et al. (I 2010al ) 
find some hints of atmospheric pollution in the photosphere of 
the host star Ibgui et al. (120101 ) suggest that the planet is even 
being disrupted by tidal forces because it may be crossing the 
Roche limit while passing the periastron. Observing WASP-12 b 
in near-ultraviolet, Fossati et al. (1201 Obi ) detect many metallic 
atoms and ions in the planetary exosphere, confirmed to exceed 
the Roche's lobe. The same conclusion is reached by Croll et al. 
( 120101) on the basis of the analysis of the near-infrared thermal 
emission of the planet. Moreover, Fossati et al. (I2010bt detect 
an early ingress and interpret this phenomenon as the presence 
of previously stripped exosphere matter, which forms a diffuse 
trail along planet's orbit. The orbital motion could produce a 
high-opacity bow shock ahead of the planet. Following this idea. 
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Vidotto et al. (1201 Oi l put an upper constraint on the strength of 
WASP-12 b's planetary magnetic field, which was found to be 
2.4 X IQ-^ T. 

Lopez-Morales et al. ([2010) have detected thermal emis- 
sion from WASP-12 b by observing its occultations in the z' 
band and support a non-circular orbit. However, Campo et al. 
(1201 0) analysed occultations in the infrared and find that they 
are centred at phase ~0.5. A model with the orbital precession 
has been proposed to explain possible variation in timing of oc- 
cultations. Then, Croll et al. (2010) observed occultations in the 
near-infrared and combined available photometric and spectro- 
scopic data to argue that the orbital eccentricity of the planet 
is indistinguishable from 0. No convincing evidence for orbital 
precession has been found. 

Li et al. (i2010 ) note that WASP-12 b's non-zero eccentric- 
ity - if real - could be excited by an additional super-Earth 
planet embedded in the circumstellar disk. This hypothesis made 
WASP-12 b an attractive target to search for transit timing vari- 
ation (TTV) caused by gravitational interactions with an inner 
perturber. 

2. Observations and data reduction 

Two transits of WASP-12 b were observed with the 2.2-m tele- 
scope at Calar Alto Observatory (Spain) during two runs on 2 
and 26 February 2010. During remaining six nights granted to 
the project, bad weather conditions did not allow us to gather sci- 
entific data. The Calar Alto Faint Object Spectrograph (CAFOS) 
in imaging mode was used as a detector It was equipped with 
the SITe CCD matrix (2048 x 2048, 24fim pixel). A 630 x 660 
pixel subframe was used during observations to shorten read-out 
time to about 30 s. The effective field of view was limited to 
5'6 X 5'8 with the scale of 0.53 arcsec per pixel, so wide enough 
to record nearby comparison stars. Including overheads, the ob- 
servations were carried out with the cadence between 77 and 95 
s, depending on the exposure time, which was refined during out- 
of-transit phase to achieve maximal efliciency. To avoid affect- 
ing transit timing, no exposure-time changes were made during 
the ingress or egress phases. The photometric monitoring was 
performed in the Johnson /?-band filter The telescope was sig- 
nificantly defocused, and stellar profiles exhibited a donut shape. 
This method was expected to minimise random and flat-fielding 
err ors (S outhwoith et al. 12009a] I2009bj I20l0i also Johnson et 
al. 120101 . The stellar flux was spread over a ring of ~24 pix- 
els (~13 arcsec) in a diameter This value allowed to avoid con- 
tamination of faint neighbour stars, the closest one of which is 
located ~9 arcsec from the target star The stellar images were 
kept exactly at the same position in the detector's matrix during 
each run thanks to auto-guiding. Our tests showed that in this 
case dividing by a flat-field frame has no detectable effect on a 
final light curve. Precise timing was assured by synchronising 
the computer's clock to Coordinated Universal Time (UTC) by 
Network Time Protocol software, accurate to better than 0.1 s. 

During the first run, sky conditions were mainly photomet- 
ric with ephemeral thin clouds before egress. A gap just be- 
fore the transit's beginning was caused by technical problems. 
In the second run, observations were gathered in gaps between 
clouds. Conditions during ingress and egress were photometric 
and partly cloudy during the flat bottom phase. Observations in 
both runs were carried out in bright time, and the-second run 
data were particularly affected by moonlight. The details on ob- 
servations are presented in Table [T] 

The CCD frames were processed using a standard procedure 
that included subtraction of a dark frame and flat-fielding with 



Table 1. The summary of observing runs. 



Run* 


Date 


N 


X 




1 


02 Feb. 2010 


167 


1.19 -> 1.01 ^ 1.16 


45, 50, 60 


2 


26 Feb. 2010 


69 


1.01 ^ 1.40 


42, 50, 60 



Notes. Ncxp ~ the number of useful exposures, X - airmass changes 
during the run, Te^p - exposure times, dates in UT at the beginning of 
nights. 



dome flats. The magnitudes of WASP-12 and comparison stars 
were determined with aperture photometry with aperture radii 
ranging from 8 to 20 pixels. The smallest photometric scatter 
was found to be at an aperture radius of 12 pixels. This value is 
also a safety limit to avoid neighbour-star contamination. As a 
result of iterative tests, two nearby comparison stars, for which 
airmass and colour differences were the smallest, were used to 
set the zero level of the photometric scale. The formal photo- 
metric error of the individual measurements was found to be 
in a range between 0.6 and 0.7 mmag. Trends in light curves, 
caused by differential atmospheric extinction, were modelled us- 
ing out-of-transit measurements and then removed using the fol- 
lowing approximation. For data from 2 February, the second- 
order polynomial was fitted with the least-square method and 
then subtracted from the light curve. The out-of-transit data on 
26 February were found to be of marginal use, so only a linear 
trend could be determined and removed. 



3. Results 

The further analysis here was inspired by South worth C2008J . 
Individual light curves were modelled with the JKTEBOP code 
(Southworth et al. 200 4a, 2004b), which is based on the EBOP 
program (Etzel 119811 Popper & Etzel 1981). Five parameters 
describing a shape of a light curve were allowed to float dur- 
ing the fitting procedure. We used fractional radii of the host 
star and planet, defined as r» = ^ and '"b = respectively, 
where and R\, are the absolute radii of the bodies and a is the 
orbital semi-major axis. In practise, the combinations of these 
parameters were used, a sum r, + rb and ratio k - r^/rb, because 
they were found to be only weakly correlated with each other 
(Southworth '2008). The directly fitted orbital inclination, /, al- 
lows calculation of the transit parameter b = cos ;. The ini- 
tial values of parameters listed above were taken from Hebb et 
al. (12009). To detect any variation in transit timing between ob- 
served transits, the mid-transit time was a free parameter whose 
initial value was calculated according to the ephemeris given by 
Hebb et al. (2009). High-quality light curves need to set limb- 
darkening coefficients (LDCs) as free parameters and test differ- 
ent limb-darkening (LD) laws (Southworth '2008). We therefore 
considered linear, logarithmic, and square-root LD, for which 
theoretical LDCs in the Johnson R band were bilinearly inter- 
polated from tables by 'Van Hamme ( 119931 ). In a first iteration, 
we allowed both linear u and non-linear v coefficients to float, 
but unrealistic values were derived. Therefore, we kept v fixed 
at theoretical values and allowed only u to vary in final fitting 
runs. The contribution of v into the error budget was included by 
perturbing it by ±0.1 around theoretical values and assuming a 
flat distribution. 

The errors of derived parameters were determined in two 
ways for each combination of a data set and adopted LD law. 
Firstly, we ran 1000 Monte Carlo (MC) simulations, and a spread 
range of a given parameter within 68.3% was taken as its error 
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Table 2. Parameters of transit light-curve modelling derived for 
light curves in runs #1 and #2. 



Parameter 



Linear 



Logarithmic 



Square-root 



Run#l 



r, + 
k 

n, 

i (deg) 
h (R,) 
u 

Tq (JDutc) 
rms (mmag) 

r, + rt, 

k 

r, 
rb 

i (deg) 
b{R.) 
u 

Tq (JDutc) 
rms (mmag) 

^red 



n I 1 qs7+0'""'53 

<J.il7JZ, „„„„„ 

n ■xAT.q+O-Oois 
04110+""°** 

"■-'-'-0.03 

0.5896 
1.1647 



^74+0.005 

"•-^ -0.005 

n 1 1 Tss+O-Oooes 

u. 1 i ' JJ_o.00072 
n T.-iAl+O-OO'H 

03928+°"°"" 

"■"-'^^°-0.p0068 

82.2+[j« 
n QQi +0.041 

571+" °^'* 

V.J I -l_o.Q54 

40146+°"°'"2 
0.5840 



1.1390 



Run#2 



379+""'^ 

n 1 I q7+?iSf)^6 
^ -0.0021 

n ^QQ+0.016 
0404+°"°^7 

"■"^"^-0.0018 
Q 4Q+Q.f4 

398+°"^*^ 
41536+°"°"" 
0.9747 
2.8718 



+0.016 



0.382 

"•ii°°-o.ooi6 
0.341+"°'" 

0406+^-^''^^^ 

81.8![| 

40+" '" 

"•^ -0.13 

57+" °^ 

—0 08 

41536+°"°"" 
0.9748 
2.8664 



0.3710+"°""'* 
1175+8:885^ 

0.3320+"°""^ 
0390+8-88^^ 

"■"-'^"-0.0007 

82.5!!J? 

n ^7^+0 043 
"•■"-'-0.039 

-0 024+" °"* 

"■"^^-0/)46 

40144+"°"°" 

"■^"^^^-0.00010 

0.5839 
1.1373 



^0^+0.017 

ii90+" °"^3 

"■'^ ^^"-0.0015 
n ^4^+0 014 

0.0407^1^^^ 
81.6il 
41+"'" 

"■^^-0.11 

-0 017+"°*" 

"■"^^'-0/146 

41536+"°"°'^ 

"■^^-^-^"-0.00014 

0.9748 
2.8660 



Common parameter 
0.2776° 



0.6931° 



Notes. To is based on UTC and is given as 2455230-t- and 2455254-)- for 
run #1 and #2, respectively. 

permuted by ±0. 1 on a flat distribution 



estimate. Secondly, the prayer-bead method (e.g. Desert et al. 
120091 Winn et al. 120091 ) was used to check whether red noise is 
present in our data sets. Errors returned by both methods were 
found to be similar for run#l, which allowed us to adopt MC 
errors as the final ones. The run#2 light curve is of lower qual- 
ity, and this was reflected in the error analysis. The prayer-bead 
method gave error values 2-3 times greater than the MC algo- 
rithm. This indicates that, besides Poisson noise, there is cor- 
related (red) noise in the light curve. In this case we adopted 
prayer-bead errors as the final ones. The only exception is the 
mid-transit time whose MC errors were found to be greater, so 
they were used as the final ones. This may suggest that the mid- 
transit time error may still be underestimated. 

Table|2]contains the results of individual fits. For lower qual- 
ity data, represented by our run#2 light curve, we found that the 
goodness of the fit does not depend on the LD law. A signifi- 
cant value of rms (0.97 mmag) clearly indicates that data are af- 
fected by non-photometric conditions during the night, so results 
should be treated with caution. Only the mid-transit time seems 
to be reliable as non-Poison noise did not aff'ect it. In the case 
of the run#l high-quality data, using the linear LD law results 
in a slightly poorer fit, while non-linear LD laws give compara- 
ble values of xl^^- The same effect was noticed for WASP-4 by 
Southworth (2009b). The light curves acquired in both runs are 
plotted in Fig.[l] together with best-fitting models and residuals. 



4. Physical properties of WASP-12 system 

Results of light-curve modelling allowed us to calculate plane- 
tary, stellar, and geometrical parameters. We used results from 
the best-fitting model based on the run#l data and square-root 
LD law. The value of was calculated assuming a - 0.0229 ± 




J I I I I I I I I I I I I I I I I I I I L 

-0.10 -0.05 0.00 0.05 0.10 

Time from mid-transit (d) 

Fig. 1. Light curves of two transits of WASP-12 b observed on 2 
and 26 Feb. 2010. The best-fitting models, based on the square- 
root limb darkening law, are plotted with continuous lines. The 
residuals are shown in the bottom panels. 



0.0008 AU (Hebb et al. |2009l ). The planetary mean density pb 
uses WASP-12 b's mass of 1.41 ±0.10 Mj, known from radial ve- 
locity measurements. The surface gravitational acceleration, g\„ 
was calculated according to a formula (Southworth et al. .2007l l: 



27r Vl -eV 



Ph 



sin / 



(1) 



where Pb is the orbital period of the planet (see Sect. 5), and 
is the stellar velocity amplitude taken from Hebb et al. (I2009I I. 
The orbital eccentricity was assumed to be 0. 

The equilibrium temperature, Teq, was derived by assum- 
ing the effective temperature of the host star Teff = 6300+^qq 
K (Hebb et al. 2009) and using the relation (Southworth ,2008i) : 



1 -A 
AF 



"■25 , >0 5 



(2) 



where A is the Bond albedo and F the heat redistribution factor. 
This formula may be simplified by assuming relation A - \ - 
AF, and the modified equilibrium temperature was derived 
(Southworth 2008). 

We calculated the Safronov number 0, which is proportional 
to the ratio of the escape velocity from the planet and velocity of 
its orbital motion. Southworth (12008) gives formula 



= 



(3) 
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Table 3. Physical properties of the WASP- 12 system derived 
from light-curve modelling (see text for explanations). 



Table 4. Results of transit timing. 



Parameter 



This work 



Hebbetal. (2009 1 





Planetary properties 




Rb (R,) 


1.9 + 0.1 


1.79 + 0.09 


Pb (Pj) 


0.22 ± 0.04 




gh (m s-2) 


10.3 + 0.6 


9.8 + 0.7 


(K) 


2570!^™ 


2516 + 36 





0.026 ± 0.005 




/ 


0.083 

Geometrical parameters 




i (deg) 


82.5!!J-? 


83.r|* 


b(R,) 


-^-0.039 


"■-'"-0.06 


(Rb/RS- 


0.01380 + 0.00016 0.0138 + 0.0002 




Stellar properties 




R. (Re) 


1.63 + 0.08 


1.57 + 0.07 


p. (po) 


4.14 + 0.09 


4.17 + 0.03 


logg. (cgs) 


0.31 + 0.04 


0.35 + 0.03 



where the the value of 1.35 + 0.14 Mq was taken as the star's 
mass M, (Hebb et al. l2009b . 

The shape of the planet was found to be far from an ideal 
sphere. The oblateness of the planet, directly returned by the 
JKTEBOP code and defined as 



1-2 



(4) 



where ri and r2 are the polar and equatorial radii, respectively, 
was found to be 0.083. This locates WASP-12 b between Jupiter 
and Saturn, whose oblateness values are 0.065 and 0.098, re- 
spectively. Results of our calculations are collected in Table [3] 
where the literature values are also given. For almost all param- 
eters we obtained excellent agreement within Icr error bars with 
determinations repotted by Hebb et al. (I2009t . 

5. Transit timing 

The mid-transit times were transformed from JD based on UTC 
into BJD based on Bary centric Dynamical Time (TDB) using the 
on-line converteiQ by Eastman et al. (2010). Having a long time 
span of observations, we redetermined an improved ephemeris. 
As a result of fitting a linear function of the epoch and orbital 
period P\,, we obtained 

To = 2454508.97682 ± 0.00020 BJDtdb 

Pb = 1.09142245 ± 0.00000033 d. 

Individual mid-transit eiTors were taken as weights. The obser- 
vation minus calculation (O - C) diagram, plotted in Fig. |2] 
clearly shows a deviation in our mid-transit times from the 
new linear ephemeris. The amplitude of the possible TTV of 
WASP-12 b was found to be 73 s between 22 epochs at the 
level of 3.4cr. To check whether keeping all parameters free for 
run#2 data may affect the mid-transit time, we repeated mod- 
elling iTin#2 data with all parameters except the mid-transit time 
kept fixed. Their values were adopted from the run#l analy- 
sis. The derived mid-transit time agreed with the full-modelling 
value well within Icr. As the additional test, both light curves 
were reanalysed with a model-fitting algorithm available via the 
Exoplanet Transit Database^ (Poddany et al. l2010l l to obtain 



Run# 


To (BJDtdb) 
Tq (JDutc) 


Tq error (d) 


Epoch 


- C (d) 


1 


2455230.40673 


+0.00011 


661 


-0.00033 




2455230.40144 








2 


2455254.41887 


+0.00014 


683 


+0.00051 




2455254.41536 









Notes. The mid-transit times, Tq, are given as both BJD (based on 
Barycentric Dynamical Time, TDB) and JD (based on Coordinated 
Universal Time, UTC). 



3 
U 

6 




-0.0006 



Fig. 2. The observation minus calculation (O - C) diagram for 
WASP-12 b generated for a new linear ephemeiis. The open 
symbol comes from Hebb et al. ( I2009I I, and the filled symbols 
denote mid-transit times reported in this paper 



mid-transit times with an alternative method. The procedure em- 
ploys the OCCULTSMALL routinc of Mandel & Agol (2002) and the 
Levenberg-Marquardt non-linear least-square fitting algorithm, 
which also provides uncertainties. The procedure assumes the 
linear LD law with the fixed theoretical LDC whose value was 
taken from tables by Van Hamme (1993). We derived To equal to 
2455230.40143 + 0.00012 and 2455254.41533 ± 0.00020 JDutc 
for tuns#l and #2, respectively. Mid-transit times obtained with 
both codes turned out to be consistent with each other within 1 cr, 
which strengthens the reliability of the possible TTV detection. 
The transit-timing results are summarised in Table 

The short time scale of possible TTV signal suggests that it 
may be caused by an additional low-mass planet perturbing the 
orbita l motion of WASP-1 2 b (M iralda-Esc udei2002 , Schneider 
120041 Holman & Murray |20E1 Agol et al. l2005l Steffen et al. 
'2007). Because of the few points in the O - C diagram, it is 
impossible to characterise the time scale of the TTV and de- 
tect some peiiodicity. However, by knowing the amplitude of the 
TTV, which may be treated as a minimal value, one can deter- 
mine the minimal mass of the additional planet in a wide range 
of orbital configurations. We generated 745 synthetic (9 - C di- 
agrams for orbits from the inner 1:5 mean motion resonance 
(MMR) to the outer 5:1 one. For simplicity we assumed that the 
planetary orbits are coplanar and initially circular. We used the 



' http://astroutils.astronomy.ohio-state.edu/time/utc2bjd.html 
^ http://var2.astro.cz/ETD/ 



Mercury package (Chambers |1999|l and the Bulirsch-Stoer al- 
gorithm to integrate the equations of motion for this three-body 
problem. Computations covered 700 periods of the transiting 
planet, i.e. the time span of observations. Figure |3]illustrates the 
results of our simulations. The observed amplitude of the TTV 
could be caused by a body as small as an Earth-mass planet if 
it was in one of several low-order MMRs. For instance, such 
a planet would induce stellar radial-velocity variations with the 
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1:2 3:5 5:3 2:1 3:1 




0.2 0.5 1.0 2.0 5.0 



Fig. 3. The mass of a perturbing planet that could generate a TTV 
signal with an observed amplitude of 73 s, as a function of ra- 
tio of orbital periods of transiting planet, P^, and the perturber, 
Pp. Vertical dashed lines indicate MMRs. Orbits in the range of 
Pp/Pb between 0.69 and 1 .45, marked by a grey area, were found 
to be unstable. The case of Pp/f b = 1 was found to be stable but 
was skipped in the plot. 

semi-amplitude of 0.6 and 0.4 m s ' in the inner and outer 2:1 
resonances, respectively, which in practice are undetectable with 
modern state-of-the-art instruments. This finding makes WASP- 
12 b an attractive target for the TTV follow-ups. 
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6. Concluding remarks 

Our high-precision photometry gathered for two transits of 
WASP- 12 b in February 2010 allowed us to redetermine param- 
eters of this extra-solar planetary system. They were found to 
agree with the values reported by Hebb et al. (2009 1. Two new 
mid-transit times were used to derive improved ephemeris. We 
found that the observed transit timing suggests there is a short- 
time variation at the level of 3 Act. Such a TTV signal could be 
caused by an additional teiTestrial-type planet if both planets are 
close to orbital resonances. Further high-precision observations 
of WASP-12 b's transits are needed to verify the existence of the 
TTV signal and put constraints on the second planet's parame- 
ters. 
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